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We propose a simple model with the 1^ symmetry in order to answer whether the symmetry 
is a good concept in QCD with light quark mass. The model is constructed by imposing the 
flavor-dependent twisted boundary condition (TBC) on the three-flavor Polyakov-loop extended 
Nambu-Jona-Lasinio model. In the model, the Zn symmetry is preserved below some temperature 
T c , but spontaneously broken above T c . Dynamics of the simple model is similar to that of the 
original PNJL model without the TBC, indicating that the In symmetry is a good concept. We 
also investigate the interplay between the symmetry and the emergence of the quarkyonic 
phase. 
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1. Introduction 

In the limit of zero current quark mass, the chiral condensate is an exact order parameter for 
the chiral restoration. In the limit of infinite current quark mass, on the contrary, the Polyakov 
loop becomes an exact order parameter for the deconfinement transition, since the symmetry is 
exact there. For the real world in which u and d quarks have small current quark masses, the chiral 
condensate is considered to be a good order parameter, but it is not clear whether the Polyakov loop 
is a good order parameter. In this paper, we approach this problem by proposing a simple model 
with the Zjv symmetry. This paper is based on our recent papers 01. 

We start with the SU(3) gauge theory with three degenerate flavor quarks. The partition func- 
tion Z in Euclidean spacetime is 



DqDqDA exp[-S ] (1-1) 
with the action 

So = j d 4 xl£q f (Y v D v + m f )q f +^FZ v 2 ] (1.2) 
and the temporal boundary condition 

= i/r) = -?/(*,()). (i.3) 

The Z3 transformation changes the ferimon boundary condition as ||3|, |j] 

q f (x,P) = -&xp(i2xk/3)q f (x,0) (1.4) 



for integers k = 0, 1,2, while the action So keeps the original form (1.2) since the Z3 symmetry is 
the center symmetry of the gauge symmetry The Z3 symmetry thus breaks down through the 
fermion boundary condition. Now we assume the twisted boundary conditions (TBC) 

q f (x,p) = -expH0 / ) 9/ (x,O) = -expH(0 1 +2^(/-l)/3)] 9/ (^O) (1.5) 

for flavors / = 1,2,3, where 6\ is an arbitrary real number in the range < 6\ < 2n. QCD with 
the TBC has the Z3 symmetry. Actually the Z3 transformation changes / into f — k, but f — k can 
be relabeled by / since So is invariant under the relabeling. The TBC is useful to understand the 
color confinement. 

When the fermion field qj is transformed as 

qf = exp (—i9fTT)q' f (1.6) 

for Euclidean time T, the action So is changed into 

5(0/) = / d 4 x[^q' f (YvD v - HfY4 + mfWf + A*?vl ( L7 > 
J f *g 

with the imaginary quark number chemical potential \if = iTdf, while the TBC is transformed back 



to the standard one (1.3). The action So with the TBC is thus equivalent to the action S(0/) with 
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the standard one (|L3|). The partition function Z(T, 6) has the Roberge-Weiss (RW) periodicity [Qj: 
Z(T,6) = Z(T,6 + 27ik/3) for any integer k. The Polyakov-loop extended Nambu-Jona-Lasinio 
(PNJL) model is a good model to understand QCD at finite imaginary chemical potential 0. 

In this paper, we consider QCD with the TBC in order to answer whether the Za? symmetry 
is a good concept in QCD with light quark mass. Dynamics of the theory is studied concretely by 
imposing the TBC on the PNJL model, i.e., the TBC model. 

2. PNJL and TBC models 

The three-flavor PNJL Lagrangian is 



& = L?/(yvOv-M/74 + m/)?/-G s ££[(?/A,^ / ) 2 + ( g/ jy 5 ^ / ) 2 ] 

/ / a=0 



-Gd 



detft(l + 75 ky + detft(l - y 5 )qj + & (4»[A] , <J>* [A] , T). 

'j 'j J 



(2.1) 



For the Polyakov potential % as a function of the Polyakov-loop <I> and its conjugate <J>*, we take 



the potential of Ref. [g]. Now we consider the PNJL model with the TBC (|L5j), that is, the TBC 
model. The thermodynamic potential of the TBC model is nothing but that of the PNJL model 
with the flavor-dependent imaginary chemical potential \ij = id/T. In the mean-field level, the 
thermodynamic potential £2 of the TBC model is obtained as 



Q 



2 E E /(Sn 

f=u,d,sc=r,g,b J V^""7 L 

+ c/(cT M ,cr d ,<T i )+^(a) ) a)*,r), 



£■/ + ^-In [1 + j*°tPf e - pE J] + ^ In [1 + 



(2.2) 



where oy = (q/q/), Ef = y P 3 +-^/ 2 with 

Mf = mf-4GsOf + 2GnOfOf» (2.3) 
for / / /' and / / /" and /" / /'". The mesonic potential [/ (a„, a^, a 4 ) are obtained by 

U(a u ,a d ,o s )= £ 2G s a)-4G D a 1( a rf a v . (2.4) 

f=u,d,s 



The vacuum term in (2.2) is regularized with the three-dimensional cutoff A. For the parameter set 
(Gs, Gd, mi, m s , A), we take the set of Ref. [g], except that the s-quark mass m s is taken to be the 
same as the light quark mass mi = m u = m d . 



Taking the color summation in (2.2), we can get 



a = -2 Y 



d ' P [ N c E f + Iln[l+C3,i(p)^+C3, 2 (p)e 2 ' e /+C 3 ,3(p)^1 



+ iln[l +CtMe- id fe-^ E f + C 3 * 2 (pK 2ie ' + C 3 * 3 (pK 
+ U(a u ,(J d ,a s ) + ^(^*,T) 



3/e, 



(2.5) 
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with 

C 3) i(p)=3*e-^, C 3 , 2 (p) = 30*e- 2 ^, C 3)3 (p) = e~^ E f. (2.6) 

When <I> = 0, £1 has no flavor dependence, since C 3) i = C3 2 = and the factors e ±3,e f do not 
depend on flavor. The flavor symmetry is thus preserved in the confinement phase with <J> = 0. 



3. Numerical results 

Figure |l| shows T dependence of order parameters a, <I> and a = o u — o d = o u — a s in (a) the 
PNJL model with {d u ,d d ,d s ) = (0,0,-0) = (0,0,0) and (b) the TBC model with (6 u ,6 d ,0 s ) = 
(0,6,— 6) = (0,2n/3, — 2n/3); note that ao is an order parameter of the flavor symmetry. In the 
PNJL model, both the chiral restoration and the deconfinement transition are crossover. In the TBC 
model, the first-order deconfinement transition takes place at T = T c w 195 MeV. Below T c , ao and 
<I> are zero, as expected. The flavor symmetry is thus preserved by the color confinement. Above 
T c , ao and <I> become finite, indicating that the flavor and Z3 symmetries break simultaneously. The 
chiral restoration is very slow in the TBC model because of the breaking of flavor symmetry. As 
mentioned above, the deconfinement transition is first-order at 6 = 2n/3, but crossover at 6 = 0. 
This means that there appears a crtical endpoint at some 6 when 6 is varied from 27i/3 to 0. 
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Figure 1: T dependence of a, 4> and a Q for (a) (8 U , d , 6 S ) = (0,0,0) and (b) (0„, d , S ) = (0,27r/3,47r/3). 
<7 is normalized by the value Oo at T = 0. See legends for the definition of lines. 

Next we consider the entanglement-PNJL (EPNJL) model gj. The four-quark vertex G§ 
is originated in a gluon exchange between quarks and its higher-order diagrams. If the gluon field 
A v has a vacuum expectation value (Ao), A v is coupled to (Ao) and hence to [§]. This effect 
allows Gs to depend on <!>: Gs = G$(^) [§]. In this paper, we simply assume the following G$(^>) 



by respecting the chiral symmetry, the charge-conjugation symmetry [|10J and the extended Z3 
symmetry ph: 



G s (4>) = G s [l - ai^O* - a 2 (<I> j + <I>* 3 )]. 



(3.1) 



In principle, Gd can also depend on <£. However, ^-dependence of Gd yields qualitatively the 
same effect on the phase diagram as that of Gs [S]- We can then neglect ^-dependence of Gd- 



The parameters ai and a 2 in (3.1) are so determined as to reproduce two results of LQCD at finite 
T; one is the result of 2+1 flavor LQCD at /I =0 [11] that the chiral transition is crossover at 
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the physical point and another is the result of degenerate three-flavor LQCD at 6 = % [12] that 
the order of the RW endpoint is first-order for small and large quark masses but second-order for 
intermediate quark masses. The parameter set {(X\ , a?) satisfying these conditions is located in the 
triangle region [H 

{-1.5ai +0.3 < a 2 < -0.86c*i +0.32, a 2 > 0}. (3.2) 

As a typical example, we take a\ = 0.25 and a 2 = 0.1, following Ref. 

In Fig. ||, a, a and <I> are calculated as a function of T with the EPNJL model for (a) 
(9u,9d,0 s ) = (0,0,0) and (b) (d u ,6 d ,d s ) = (0,27r/3,47r/3). In panel (a), the chiral restoration 
and the deconfmement transition are first-order, because the current quark mass (5.5MeV) is small 
and the correlation between ay and <I> is strong [@]. In panel (b), one can see that <J> = a = in the 
confinement phase. The EPNJL model with the TBC yields similar T dependence to that without 
the TBC for both the chiral restoration and the deconfmement transition, since the flavor-symmetry 
breaking above T c is weakened by the strong correlation between ay and <£. 
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Figure 2: T dependence of order parameters a, ao and <t> calculated with the EPNJL model for (a) 
(0 u ,G d ,9 s ) = (0,0,0) and (b) (9 u ,G d ,9 s ) = (0,2^/3,4^/3). Here (J and a are normalized by a . Note 
that «o = in panel (a) and ao > in panel (b), while a < in both panels. See legends for the definition of 
lines. 

Now we consider the PNJL model with the flavor-dependent complex chemical potentials 
iif = il + iTdf in which the 0/ are defined by 



(0/) = (0,0,-0) 



(3.3) 



The present model with the \if is reduced to the standard PNJL model with the flavor-independent 
real chemical potential \i when 0=0 and to the TBC model with real \i when = 27T/3. Varying 
from Oto0=27r/3, one can see how the phase diagram changes between the approximate 
color-confinement in the standard PNJL model and the exact color-confinement in the TBC model. 

Figure || shows the phase diagram in the \i-T plane. Panels (a)-(c) correspond to three cases of 
0=0, 87T/15 and 2n/3, respectively. When 0=0, both the chiral and deconfmement transitions 
are crossover at smaller }i, but the chiral transition becomes first-order at larger }X. When = 27r/3, 
the deconfmement transition is the first-order at any ju, whereas the chiral transition line becomes 
first-order only at H w M f = 323 MeV. The region labeled by "Qy" at /I > M f and small T is the 
quarkyonic phase [13], since = but the quark number density n is finite there. The region 
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labeled by "Had" is the hadron phase, because the chiral symmetry is broken there and thereby the 
equation of state is dominated by the pion gas The region labeled by "QGP" corresponds to 
the quark gluon plasma (QGP) phase, although the flavor symmetry is broken there by the TBC. As 
6 decreases from 2n/3 to zero, the first-order chiral transition line declines toward smaller /I and 
the critical endpoint moves to smaller jU. Once 6 varies from 2n/3, the quarkyonic phase defined 
by <I> = and shrinks on a line with T = and jj. > Mf and a region at small T and /I > Mf 
becomes a quarkyonic-like phase with small but finite <I> and n ^ 0; the latter region is labeled by 
"Qy-like". 
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Figure 3: Phase diagram in the ji-T plane. Panels (a)-(c) correspond to three cases of 9 = 0, 87T/15 and 
2n/3, respectively. The thick (thin) solid curve means the first-order deconfinement (chiral) phase transition 
line, while the thick (thin) dashed curve does the deconfinement (chiral) crossover line. The closed circles 
stand for the endpoints of the first-order deconfinement and chiral phase transition lines. In panels (a) and 
(b), the thick-solid line at T = and > Mf = 323 MeV represents the quarkyonic phase. 



4. Summary 

We have proposed a simple model with the Zjy symmetry in order to answer whether the 
symmetry is a good concept in QCD with light quark mass. The model called the TBC model is 
constructed by imposing the flavor-dependent twisted boundary condition ( |L~5| ) on the PNJL model. 

In the TBC model, the Z3 symmetry is preserved below T c , but spontaneously broken above 
T c . Below T c , the color confinement preserves the flavor symmetry. Above T c , meanwhile, the 
flavor symmetry is broken explicitly by the TBC. The flavor-symmetry breaking makes the chiral 
restoration slower, but the entanglement interaction between a and <J> makes the restoration faster. 
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Consequently, we can except that QCD with the TBC is similar to original QCD with the standard 
quark boundary condition. 

We have also investigated the inteiplay between the symmetry and the emergence of 
the quarkyonic phase, considering the complex chemical potentials \if = \i + iT6f with (6/) = 
(0,6,-d) in the PNJL model. The PNJL model with the \x f is reduced to the PNJL model with 
real /I for 6 = 0, but to the TBC model with real /I for 6 = 2n/3. When 6 = 2%j3, the quarkyonic 
phase defined by <I> = and n > really exists at small T and large fi. Once 6 varies from 2n/3 to 
zero, the Zn c symmetry is broken and thereby the quarkyonic phase exists only on a line of T = 
and /i > Mf. The region at small T and large /I is dominated by the quarkyonic-like phase char- 
acterized by small but finite <I> and n > 0. The quarkyonic-like phase at 6 = is thus a remnant 
of the quarkyonic phase at 6 = 2n/3. Since the Zjv symmetry is explicitly broken at 6 = 0, it is 
natural to expand the concept of the quarkyonic phase and redefine it by a phase with small <J> and 
finite n. For this reason, the quarkyonic-like phase is often called the quarkyonic phase. The gross 
structure of the phase diagram thus has no qualitative difference between 6 = 2n/3 and zero, if the 
concept of the quarkyonic phase is properly expanded. In this sense, the Z3 symmetry is a good 
approximate concept for the case of 6 = 0, even if the current quark mass is small. 
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